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. Abstract. Weak lensing observations and supernova observations, combined 

j^jy with CMB observations, can both provide powerful constraints on dark energy 

, properties. Considering statistical errors only, we find luminosity distances in- 

■ ferred from 2000 supernovae and large-scale (I < 1000) angular power spectra in- 

ferred from redshift-binned cosmic shear maps place complementary constraints 
, on wo and w a where w{z) — wq + w a (a — 1). Further, each set of observa- 

tions can constrain higher-dimensional parameterizations of w(z); we consider 
eigenmodes of the w{z) error covariance matrix and find such datasets can each 
constrain the amplitude of about 5 w(z) eigenmodes. We also consider another 
parameterization of the dark energy. 



1. Introduction 



In the past half-decade a great variety of methods have emerged for 'observing 
dark energy', enough to fill an entire 4-day meeting. One might ask why we 
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• Independent methods provide the ultimate systematic error test. 



need such a great variety of methods. To this question we give three answers: 



• Independent methods provide for less model-dependent probes of cosmic 
^ , acceleration. 

• Different methods constrain different directions in the parameter space. 

Given the importance of the dark energy mystery and the challenges to 
constraining its properties, this diversity of methods is a blessing. In this talk 
we concentrate on two methods: weak lensing and supernovae. Both of these 
highly different types of observations are potentially powerful probes of dark 
energy. Here we examine item 3: the complementary nature of the statistical 
errors. 

Weak lensing observations can deliver enormously rich datasets; there are a 
multitude of ways of using these data to constrain dark energy. Here we concen- 
trate solely on what can be done with the data on large scales because this is sim- 
plest to model, most robust to increases in shape noise and spurious psf power, 
and because others at this meeting (Bernstein and Jain) discussed the potentiall y 
powerful uses of sma l ler-scale data. See , for example. iRefregier et al. 
iTakada fc .Tainl (EiMslblTvson et al.1 (I2003T ) : I.Tain fc Tavlod (feOO.Sl h lBernstein fc .Tainl 
(|20r)4h : ISong fc Knoxl (|2003h : lHu fc .Tainl ^200.^ . 
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In Section 2 we present our models of the datasets and in Section 3 the 
constraints on wq and w a . In Section 4 we examine constraints on higher- 
dimensional parameterizations of w{z) as well as discuss another parameteriza- 
tion of the dark energy. In Section 5 we conclude. 



2. Models of the Data 

In this talk we assume a CMB survey, a cosmic shear survey and a supernova 
survey are used to constrain a large parameter space. Although we are only 
interested in the dark energy parameters here, the cosmic shear survey in par- 
ticular is sensitive to many other parameters and so we must simultaneously fit 
for them as well. The CMB survey is very useful for constraining these n on-dark- 
energy parameters. We assume Planck (comprehe nsively describ ed in ((Tauberl 
1200 j l and, more specifically for our calculations, in ISong fc Kno^ (J2003)) as our 
CMB survey since we expect these dat a to be significantly more constraining 
than WMAP data (|Bennett et al.ll2003h (due mostly to Planck's higher angular 



resolution) and these data will be available by the time we have any cosmic shear 
and super nova datasets like those we d escribe below. 

As in lKaplinghat et al.1 (I2003T ) and lSong fc Knoxl {2003) we take our (non- 
w{z)) set to be V = {uj m ,uji ) ,uj u ,9 s , z r {,k' i P^(kf),n s ,n' s ,yue}, with the assump- 
tion of a flat universe. The first three of these are the densities today (in 
units of 1.88 x 10 _29 g/cm 3 ) of cold dark matter plus baryons, baryons and 
massive neutrinos. We assume two massless species and one massive species. 
The next is the angular size subtended by the sound horizon on the last- 
scattering surface. The Thompson scattering optical depth for CMB photons, 
r, is parameterized by the redshift of reionization z x \. The primordial poten- 
tial power spectrum is assumed to be k 3 P$(k) = k^P^{kf){k/kf) ns ^ l+n 's ln ^ k / k ^ 

with kj = 0.05 Mpc -1 . The fraction of baryonic mass in Helium (which af- 
fects the number density of electrons) is yu e - We Taylor expand about V = 
{0.146,0.021,0,0.6,6.3,6.4 x 10 -11 , 1, 0, 0.24}. The Hubble constant for this 
model is h = 0.655 where H$ = 100/ikmsec -1 Mpc -1 . 



2.1. Cosmic Shear Data Model 

We call our fiducial weak lensing survey G2-7r because we imagine it as a ground- 
based survey of half of the sky. We assume a galaxy redshift distribution for a 
limiti ng magnitude in R of 26 inferred from observations with the Subaru tele- 
scope |Nagashimj^^d] ((2002). The shape of this distribution is well-described 
by the following analytic form: 



dn/dz oc z 1 ' 3 exp — (z/1.2) 



,1.2 



for z < 1 



dn/dz oc z i ' i exp — (z/1.2) 1 ' 2 for z > 1. 



(1) 



We use this distribution with the modification that half of the galaxies in the 
1.2 < z < 2.5 range are discarded as undetectable. The amplitude of the 
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distribution is such that, after this cut, the number density of galaxies is 65 per 
sq. arcmin 1 . 

We further assume that the galaxies can be divided, by photometric redshift 
estimation, into eight different redshift bins: [0-0.4], [2.8-3.2] and that for 
40 < I < 1000 systematic errors are negligible. Note that this last assumption 
is a very strong one and much work will be necessary to make it a valid one. 
Finally, we assume that the shape noise (expressed as a per-component rms 
shear) is given by 7 rms (z) = 0.15 + 0.035z. 

For more details of the data modeling, see lSong Knoxl ((20031 ) . 

The eight shear auto power spectra can be determined with high accuracy 
over a large range in I, as is shown in Figure 12.1.1 In addition to the auto power 
spectra shown, there are also 9(9-l)/2=36 cross spectra that one can measure. 
We include the cross spectra in our parameter error forecasts. They do not add 
much statistical weight because their statistical errors are highly correlated with 
the errors in the auto spectra. The large number of largely redundant 2-point 
functions will be useful f or revealing the contam inating influence of systematic 
errors. See, for example. iTakada fe White! (|2003h . 
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Figure 1. The shear-shear auto power spectra. The 8 solid curves are the 
shear power spectra from each of the galaxy source planes, B\ to B$. Dotted 
curves are the linear perturbation theory approximation. From bottom to 
top the source plane redshift ranges are By. z G [0.0,0.4], B^: z G [0.4,0.8], 
B 3 : z G [0.8,1.2], By. z G [1.2,1.6], B b : z G [1.6,2.0], B 6 : z G [2.0,2.4], B 7 : 
z G [2.4,2.8] and J5g: z G [2.8,3.2]. The error boxes are forecasts for G47T 
(see Table I) . The top dashed curve is the shear power spectrum for the CMB 
source plane. The error boxes are forecasts for CMBpol. 



1 J.A. Tyson, private communication 
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2.2. Supernova Data Model 

For the supernov a survey we assume 2000 distributed in redshift as described in 



.tor tne supernova survey we assume zuuu distributed m redsmit as described m 
iKim et al.l (|2004ft baseline SNAP supernova survey. In addition, we assume 



measurement of 100 local supernovae. To our parameter set, detailed above, 
we add a supernova luminosity calibration parameter. We do not explicitly 
include systematic errors in our analysis. The supernova forecasts in some sense 
do include systematic error estimates, because the SNAP baseline survey, by 
design, does not attempt to reduce statistical errors far below the expected 
systematic error limits (by intentionally restricting the number of supernovae in 
each redshift interval). 



3. Constraints on wq and w a 

For dark energy models with sound speeds near unity, the Jeans scale is on the 
horizon and we can ignore clustering of the dark energy on scales smaller than 
the horizon, and certainly at I < 40, the lowest I value we consider. Thus the 
dark energy results in a time-dependent but scale-independent suppression of 
power. However, the scale-independence in three dimensions does not translate 
into a scale-independence of the two-dimensional projection to shear power spec- 
tra because the projection folds the time and scale-dependence together. Note 
though that in the special case of a power-law power spectrum the projection 
does preserve the scale-dependence. Thus, in the left panel of Figure ED we see 
changing w changes the shape of the power spectra, but not at small scales 
where the three-dimensional power spectrum is well-approximated as a power 
law. In addition to altering the growth rate, dark energy alters the shear power 
spectra by changing Da{z) and therefore the projection of the fluctuation power 
at different redshifts onto the sky. 

The trend with redshift in the left panel of Figure ED can be understood as 
arising from the dependence of the growth rate on wq. Recall that we are chang- 
ing wq while keeping the angular size of the sound horizon on the last-scattering 
surface (and hence the angular-diameter distance to last-scattering) fixed. In- 
creasing wq from wq = — 1 means more dark energy at high redshift (since the 
dark energy density is now increasing with redshift as opposed to being constant) 
and less dark energy at low redshift (in order to keep the angular-diameter dis- 
tance fixed). Thus, growth is suppressed at high redshift and enhanced at low 
redshift. 

The scale-dependence of the shear power spectra dependence on the dark 
energy is a boon. The scale-dependence prevents the dark energy effects from 
being completely degenerate with a redshift-dependent calibration error. We 
have ignored shear cali bration unce r tainty in our analysis. It would be interest- 
ing to include it (as in llshak et alJ (2004)) and see how effectively the dataset 



can simultaneously constrain calibration and cosmological parameters, just with 
the assumption that the calibration is scale- independent. 
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Figure 2. Left panel shows difference between w — — .9 and w = — 1 shear 
power spectra for source bins centered at z = 0.2, 0.6, 1.4 and 3.0 from top to 
bottom. The z = 0.2 curve has no error bars on it because their large extent 
would clutter the graph. The right panel shows one a error contours in the 
wo-Wa plane for G27r, 2000 SNe and the combination (as labeled). The dashed 
curve is for G2tt with the source density uniformly decreased by a factor of 2. 

The supernova data constrain the dark energy through the dependence of 
the luminosity distance (equivalent to Da{z) in a flat Universe) on the dark 
energy. The resulting constraint from the 2000 supernovae are also shown in 
Fig. |23 We see that the two error ellipses are not perfectly aligned. Combining 
the datasets results in a factor of 2 decrease in the w a direction and a factor of 
3 decrease in the wq direction. 

An advantage of using only the large angular scale statistics is their robust- 
ness to increases in shape noise. At I < 1000, for most of the source redshift 
bins, the G27T shear errors are dominated by sample variance not shape noise 
variance. Thus an increase in the shape noise variance by a factor of 2 for all 
source bins (as would happen if the source density were uniformly lower than 
expected by a factor of two) only leads to a small increase in the error contours, 
as shown with the dashed contours in Figure 

4. Constraints on higher-dimensional parameterizations 

To deepen our understanding of how these surveys are constraining dark energy, 
we have examined how they constrain the function w(z), rather than its simple 
parameterization by wo and w a . We proceed by binning w(z) in redshift bins 
and then identifying the eigenmodes and eigenv alues of the binned w ( z) err or 
covariance matrix as was done for supernovae bv lHuterer fc Starkmanl ()200.^ . 

In the left panels of Figure we plot the three eigenmodes with the best 
determined eigenvalues. We see a striking difference in the modes for 2000 SNe 
vs. the modes for G27t: those for G2-7T stretch out to higher redshifts. The 
reason for this is that lensing is less sensitive to the growth factor at the lower 
redshifts where the source density in a given redshift bin is small and the lensing 
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Eigenvalues and Eigenmodes of w(z) 
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Figure 3. Eigenvalues (right) and first three eigenmodes (left) of the w(z) 
error covariance matrix for G27r+Planck and 2000 SNe+Planck. The large 
contributions to the eigenmodes from the highest redshift-bin are an artifact 
of that bin being much broader than the rest, extending all the way to the 
last-scattering surface. 
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window (for sources at higher redshift) is also small. Thus the supernovae are 
better at detecting changes in w(z) at lower redshift and G2-7T tends to be better 
at detecting changes at higher redshift. 

G27T and 2000 SNe also have strikingly different eigenvalue spectra. The 
error on the amplitude of the best determined mode is quite similar for each 
(~ 0.03). But the 2000 SNe spectrum is much steeper. G27T has seven modes 
with errors smaller than unity, whereas 2000 SNe has four. 

To further explore the power of these experiments to constrain dark energy 
we chose eigenmodes of w(z) because it was calculationally convenient for us. 
There are probably parameterizations with superior qualities. One we have be- 
gun exploring is modes of Pd. e .( z )> with the constraint that one of the modes 
is constant; i.e., independent of z. This approach has two virtues. First, the 
density is more directly related to the data than the equation-of-state param- 
eter. Second, interpretation of data as evidence for non-cosmological-constant 
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behavior is straightforward; one only need assess how significantly non-zero the 
amplitudes of the non-constant modes are. 

5. Conclusions 

As we have seen from Bernstein's talk and Jain's talk, weak lensing datasets 
are very rich, with dark energy constraints possible from a variety of statistical 
measures. Here we have studied just one statistic, the shear-shear two-point 
functions on large angular scales. Even restricting ourselves in this manner, 
we find the potential power of large-area weak lensing surveys is extraordinary. 
Further, we see that the statistical weight of the data constrain directions in 
parameter space different from those constrained by supernovae; i.e., they are 
complementary. Therefore the combination is particularly powerful. 

To more fully understand how these different observations probe dark en- 
ergy we examined eigenmodes of the w(z) error covariance matrix. We saw that 
the G27r observations are better at probing w{z) at high redshift and the super- 
nova observations are better at lower redshifts. A handful of eigenmodes could 
be determined with less than 100% errors. 

This work has illustrated the different statistical uncertainties on dark en- 
ergy parameters for two types of probes. A side by side comparison of compre- 
hensive error forecasts is a much more ambitious undertaking that would require 
detailed consideration of systematic errors and the extent to which they can be 
controlled. 
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